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摘 要 : 探究 昆仑 山 冰 湖 变化 特征 对 区 域 生 态 环境 和 发 展 至 关 重 要 。 基于 Google Earth Engine (GEE ERSE A 


进行 监督 分 类 ,研究 昆仑 山 近 20 a 的 冰 湖 分 布 面积 和 数目 


析 。 结 果 表明 : 
(2 


WV 


LAR ALL ,并 结合 气温 、 降 水 和 冰川 面积 进行 驱动 力 分 
1) 2000 一 2020 年 昆仑 山地 区 冰 湖 数量 增加 39.25% , 面积 增加 81.35% , 呈 西 多 东 少 的 分 布 特征 。 
昆仑 山地 区 面积 小 于 0.1 km 的 冰 湖 对 气候 变化 更 敏感 ,增长 速度 最 快 ;昆仑 山 冰 湖 主要 集中 在 海拔 4600~5600 


m, 冰 湖 数量 和 面积 分 别 占 总 量 的 71.58% 和 70.51%。(3) 2000 一 2020 年 昆仑 山地 区 气温 降低 3.45% ,降水 减少 
6.27% ,冰川 面积 减少 了 21.15% ,冰川 融化 产生 的 冰川 融 水 是 冰 湖 增长 的 主要 原因 。 研 究 结果 可 为 干旱 区 水 资源 的 


保护 和 利用 、 灾 害 预 警 等 方面 提供 科学 支撑 。 


关键 词 : GER, 昆仑 山 冰 湖 ，CEE; 气候 变化 ， 驱 动力 分 析 


冰 湖 是 以 现代 冰川 融 水 为 主要 补给 源 或 在 冰 
磺 歼 洼地 内 积 水 形成 的 天 然 水 体 "。 冰 湖 对 气候 变 
化 高 度 敏感 ,是 全 球 气候 变化 监测 的 重要 组 成 部 
分 。 古 冰川 和 现代 冰川 在 亚洲 高 海拔 地 区 造就 了 众 
多 冰 湖 ,在 全 球 气候 变 暖 的 背景 下 ,中国 西北 地 区 从 
20 1l 23 80 年 代 开 始 有 了 上 暧 干 向 暧 湿 转 变 的 趋势 ， 
Vg Jb ies YAR: AC HA DX) | PS AS BH Be Af ETF A S 
势 ,可 能 导致 冰 湖 省 决 洪水 及 泥石流 等 灾害 的 发 
^E ,对 下 游人 类 活动 构成 潜在 威胁 。 因 此 ,探讨 河 
源 区 冰 湖 变化 特征 及 影响 对 于 了 解 冰冻 圈 气 候 环 
境 变化 及 灾害 风险 评估 有 着 重要 意义 。 

目前 ,中 国 开展 的 冰 湖 变化 遥感 监测 主要 集中 
于 喜马拉雅 山 " 及 其 子 区 域 ( 叶 如 藏 布 流域 ”” \ 波 
曲 河流 域 ”\ 科 西河 流域 ”\ 希 夏 邦 马 峰 “”) CREE 
唐古拉 山子 区 域 ( 念 青 唐古拉 山西 段 ”、 帕 隆 藏 布 
TE Eug) 青藏 高 原 ( 布 加 岗 日 那曲””、 尼 都 藏 
布 流域 "彼得 藏 布 流域 ”) 、 帕 米尔 高 原 ”、 喀 喇 
昆仑 山 ”、 天 山 山脉 小 阿尔 泰山 等 区 域 。 对 于 
冰 湖 人 研究 的 侧重 点 主要 包括 :基于 遥感 的 冰 湖 信息 
提取 方法 及 应 用 一 , 冰 湖 对 区 域 气候 变化 响应 的 
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时 空 特征 分 析 及 影响 因素 分 析 *“”, 典 型 区 域 潜在 
危险 性 冰 湖 的 识别 .特征 分 析 与 预测 *, 冰 湖 演 决 
灾害 风险 评估 与 预警 等 ”””。 研 究 表 明 ,喜马拉雅 
山区 新 生 冰 湖 扩 张 速率 为 每 年 0.96 km ,其 中 冰川 
补给 湖面 积 增长 占 增 长 总 量 的 75.4% ,证 实 了 冰川 
融 水 在 冰川 湖泊 扩张 中 的 主导 作用 路 。 自 20 世 纪 
80 年代 全 球 变 暧 以 来 ,青藏 高 原 地 区 冰川 一 直 处 于 
负 平 衡 状态 ,1990 一 2010 年 青藏 高 原 冰 湖面 积 增加 
128.5 km", [HR] FE SZ f| Coe RR IER PKT B AR IRE 
Te] ,天 山 冰 湖面 积 每 年 平均 以 0.689 km? BK 0.896 85] 358 
Pe, 

昆仑 山 位 于 中 国 西北 部 的 干旱 半 干 旱地 区 , 跨 
经 度 范围 大 ,气候 要 素 变化 具有 区 域 性 特征 ,开展 
昆仑 山 冰 湖 变 化 监测 一 方面 可 以 为 干旱 区 水 资源 
的 合理 开发 利用 提供 参考 , 另 一 方面 可 为 冰 湖 省 决 
灾害 风险 评 佑 提供 数据 支撑 。 当 前 对 昆仑 山 整体 
冰 湖 变化 特征 缺乏 全 面 认识 ,因此 ,本 文 基 于 GEE 
以 2000 一 2020 年 的 Landsat 遥感 数据 开展 昆仑 山 冰 
湖 遥 感 监测 ,分 析 近年 来 昆仑 山 冰 湖 时 空 变化 特 
征 , 同 时 结合 气象 数据 ,探究 冰 湖 变化 的 影响 因素 。 
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1 研究 区 概况 


昆仑 山西 部 与 帕 米 尔 高 原 相 接 ,东部 以 青海 省 
东南 部 为 界 ,北部 与 塔里木 盆地 和 柴 达 木 贫 地 相 
XE ,南部 与 西藏 自治 区 西北 部 相 接 (图 1)。 昆 仑 山 
脉 是 一 个 东西 走向 的 山脉 ,位 于 34?~40?N,75?~ 
100°E ZZ JH] , 44€ 24 2500 km, 宽 130~200 km, 相 比较 
而 言 ,西部 比 东部 窗 ,平均 海拔 超过 了 4000 m ,总 面 
fRZ^ 50x10* km5。 昆 仑 山脉 西高 东 低 ,地形 以 较 
为 陡峭 的 山坡 和 较为 高 大 的 山峰 为 特点 ,存在 较 多 
山体 与 和 斜坡。 虽然 昆仑 山 冰 湖 接纳 一 些 雨 水 ,但 其 
主要 巾 积 雪 和 冰川 供水 ,流量 遂 因 季 而 异 , 夏 季 冰 
雪 的 强烈 融化 是 冰 湖 水 补给 的 主要 来 源 。 研 究 区 
东西 跨度 大 ,气温 .降水 等 气候 条 件 差异 较 大 ,为 便 
于 观察 冰 湖 对 气候 的 响应 ,研究 常 把 昆仑 山 分 为 
ARP 西 三 部 分 。 由 于 中 部 .东部 地 区 冰 湖 分 布 较 
少 , 本 文 在 张 连 成 等 吕 的 研究 基础 上 ,将 昆仑 山中 
部 与 东部 合并 ,以 81°E 为 界线 把 昆仑 山地 区 分 为 东 
西 两 部 分 。 


2 数据 来 源 与 方法 


2.1 数据 来 源 

遥感 影像 数据 来 源 于 GEE 平 台 , 获 取 了 2000 一 
2020 年 逐年 Landsat 系列 卫星 数据 ,用 于 提取 研究 
区 的 冰 湖 轮 廊 范围。 同时 ,为 了 减 小 季节 差异 对 
冰 湖 面积 变化 带 来 的 影响 ,降低 山区 因 云 量 过 高 而 
对 遥感 影像 解 译 造成 的 误差 ,确保 图 像 的 分 辩 率 ， 
故 大 部 分 选取 云 量 较 少 、 冰 湖 变化 稳定 的 9 一 10 月 
影像 进行 处 理 H。 其 中 ,2000 年 、2004 一 2007 年 、 
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2013 年 由 于 数据 缺失 扩展 了 影像 的 采集 时 间 (8 一 
11 月 )。 

DEM 数字 高 程 数据 采用 覆盖 中 国 昆 仑 山地 区 
范围 内 的 数字 高 程 数据 产品 SRTM DEMUTM ,数据 
的 空间 分 辨 率 为 90 m, 可 以 满足 山区 复杂 地 形 条 件 
下 的 DEM 需 求 , 用 于 提取 人 研究 区 的 海拔 、 坡 向 、 山 体 
阴影 等 地 形 因 子 , 从 而 对 冰 湖 分 布 进行 空间 分 析 ， 
获取 自 中 国 地 理 空间 数据 云 (http://www.gscloud. 
cn/)。 

气象 观测 数据 选取 FLDAS 数据 集 (https://disc. 
gsfc.nasa.gov/) 2000—2020 年 近 20 a 的 气温 和 降水 
数据 ,时 间 分 辩 率 为 1 个 月 ,空间 分 辩 率 为 0.1?x 
0.1°。 

冰川 冰 湖 编目 数据 使 用 GLIM 全 球 冰川 编目 数 
据 以 及 2017 一 2018 年 中 国 西北 冰川 编目 数据 ”, 空 
间 分 辩 率 为 30 m。 冰 湖 编 目 数据 使 用 高 亚洲 冰 湖 
编 日 数据 集 (http://www.ncde.ac.cn)。 

2.2 研究 方法 

2.2.1 基于 CEE 过 感 云 计算 平台 的 影像 监督 分 
类 本 文 基 于 GEE 平 台 的 Landsat-5/7/8 遥感 影 
像 ,通过 去 云 函数 进行 无 云 影像 的 合成 ,以 保证 研 
究 区 影像 质量 。 其 次 ,在 合成 影像 上 选取 一 定数 量 
均匀 随机 分 布 的 样本 点 ,分别 为 植被 . 裸 地 、 积 雪 、 
冰 湖 4 类 。 和 筛选 其 中 70% 的 样本 ,使 用 Random For- 
est (RE ,随机 森林 ) 进 行 影像 分 类 ,另外 30% 样 本 点 
用 于 精度 验证 (图 2)。 随 机 森林 法 具有 人 工 干预 
少 分 类 效果 明显 .运算 效率 高 等 优点 。 构 建 混淆 
矩阵 ,进行 精度 检验 ,计算 总 体 分 类 精度 (Overall 
Accuracy , OA ) 平 均值 为 0.92, Kappa 系数 平均 值 为 
0.89。 将 提取 得 到 的 冰 湖 边界 进行 矢量 化 ,进行 人 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2019)3333 号 ,对 底 图 边界 无 修改 。 下 同 。 
图 1 研究 


区 示意 图 


Fig. 1 Map ofthe study area 
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Lunas 对 蜡 常 值 不 敏感 ” ,Mann-Kendall 趋势 检验 统计 量 $ 
TMVETMHOLI 影 像 、 EH 
计算 方法 为 : 
n-l n 
Engine S= ` 3 sgn(x, - xj) (2) 
i=1j=itl 
| 监督 分 类 (随机 森林 ) | 其 中 ， 
+1, x, —%;> 0 
冰川 编目 数据 对 比 sgn(x,—x,)=40, 2,770 (3) 
-1, $44 0 
ArcGIS 
选取 时 间 序 列 长 度 w=21(2000 一 2020 年 ) ,统计 
建立 10 km 缓冲 区 


边界 拓扑 检查 冰 湖 编目 数据 对 比 


遗漏 误 判 检查 


研究 区 冰 湖 数据 


图 2 冰 湖 提取 流程 图 
Fig.2 Flow chart of glaclake extraction 


工 目 视 解 译 ,以 消除 山体 阴影 造成 分 类 误差 产生 的 
影响 。 

从 基于 GEE 得 到 的 影像 监督 分 类 结果 中 提取 
出 冰川 ` 冰 湖 , 将 得 到 的 冰川 数据 与 2015 年 冰川 编 
目 数据 进行 对 比 ” ,得 到 具有 可 靠 性 的 冰川 数据 后 ， 
以 冰川 末端 为 中 心 建 立 10 km 缓冲 区 5。 对 冰 湖 数 
据 进行 目 视 解 译 ,以 初步 消除 山体 阴影 。 之 后 与 
DEM 生成 的 坡度 图 和 冰川 缓冲 区 进行 准 加 分 析 , 进 
行 二 次 精 提 取 ,并 删除 面积 小 于 0.01 km 的 错 分 冰 
湖 ”。 提 取出 的 冰 湖 与 高 亚洲 冰 湖 编目 数据 进行 
对 比 , 冰 湖面 积 误差 为 8.69% , 冰 湖 数量 误差 为 
9.07% ,在 可 接受 范围 之 内 。 
2.2.2 Sen + Mann- Kendall Theil- Sen 斜率 估算 和 
Mann- Kendall 显著 性 检验 是 两 种 非 参数 检验 方 
IE, Sen 趋势 度 是 经 过 计算 序列 的 中 值 ,其 可 以 
很 好 地 减少 噪声 的 干扰 。 计 算 公 式 为 : 


B= wein ) Vi»i (1) 


Xx. 
i 


使 用 趋势 度 8 来 判断 时 间 序 列 的 升降 。 8>0 
时 ,时 间 序 列 呈 上 升 趋势 ,反之 , 呈 下 降 趋 势 。 
而 Mann-Kendall 方 法 本 身 对 序列 分 布 无 要 求 且 


量 $ 近 似 服 从 标准 正 态 分 布 ,使 用 检验 统计 量 Z 进 
行 趋势 检验 ,检验 中 取 显 著 水 平 c=0.05，D .= 
Zus=1.96 。Z 值 由 以 下 公式 计算 : 
pM. S»0 


[n(n 1)(2n+5)- Yui- 120, 5] 


var(S) = i8 


(5) 
式 中 :n 是 序列 中 数据 个 数 ,m 是 序列 中 结 (重复 出 
现 的 数据 组 ) 的 个 数 , t, 是 结 的 宽度 (第 i 组 重复 数 
据 组 中 的 重复 数据 个 数 )。 对 Z 值 进行 双边 趋势 检 
验 。 在 Matlab 中 运用 该 方法 完成 对 研究 区 气温 和 
降水 的 趋势 分 析 和 显著 性 检验 。 
2.2.3 Person 相关 系数 Person 相关 系数 又 称 为 
Person 积 矩 相关 系数 ,用 于 度量 两 个 变量 并 和 了 之 
间 的 线性 相关 性 ,通常 用 + 来 表示 ,r 为 两 变量 间 协 
方差 和 标准 差 的 商 , 公 式 为 : 


n 


2 (X -Xyv,- Y) 
r= —— = (6) 


n 


Se _x) pac -py 

Person 相关 系数 进行 相关 性 分 析 的 约束 条 件 为 
两 变量 独立 ;两 变量 均 为 连续 变量 ;两 变量 均 符合 
正 态 分 布 且 二 元 分 布 也 符合 正 态 分 布 ;两 变量 之 间 
有 线性 关系 。 + 什 为 [-1,1] 区 间 , 其 绝对 值 越 大 , 变 
量 之 间 的 相关 性 就 越 强 ;1=0 时 ,表示 两 变量 不 存在 
线性 相关 ;1=1 时 ,表示 两 变量 呈 完 全 正 线性 相关 ; 
,=-1 时 ,表示 两 变量 旦 完全 负 线 性 相关 。 
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3 结果 与 分 析 


3.1 昆仑 山 冰 湖 变化 时 空 特征 分 析 

3.1.1 昆仑 山 冰 湖 数量 和 面积 总 体 变 化 特征 WT 
20 a 昆 仑 山地 区 冰 湖 总 体 表现 出 数量 增多 .面积 增 
大 的 趋势 , 且 二 者 的 变化 表现 出 一 致 性 (图 3)。 
2000 年 该 地 区 共有 冰 湖 642 个 , 冰 湖 面积 共计 
28.03 km? ;2020 年 冰 湖 数量 增长 至 894 个 ,总 面积 达 
50.83 km。 冰 湖 的 总 数量 增加 了 39.25% ,总 面积 增 
加 了 81.35% ,平均 以 每 年 2.22% 和 3.65% 的 速度 增 
加 。 昆 仑 山 冰 湖面 积 平均 每 年 扩张 1.09 km ,与 高 
亚洲 冰 湖 变化 趋势 一 致 ,扩张 速率 低 于 帕 米 尔 高 原 
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积 年 均 扩张 速率 (6.43 km-a") 2 。 

昆仑 山 冰 湖 的 变化 整体 可 分 为 两 个 阶段 ,第 一 
个 阶段 是 2013 年 之 前 , 冰 湖 总 量 较 少 、 总 面积 较 小 ， 
数量 年 际 变化 明显 ,出现 周 期 性 的 先 升 高 后 降低 现 
象 。 如 ,2000 一 2004 年 .2004 一 2009 年 .2009 一 2013 
年 , 冰 湖 数量 增加 ,后 经 历 2~3 a 的 回落 ,又 降低 至 
升 高 前 水 平 。 第 二 个 阶段 是 2014 一 2020 年 , 冰 湖 数 
量 在 2014 年 大 幅 增 加 后 持续 稳定 增长 。 冰 湖面 积 
无 显著 变化 ,是 平稳 增加 趋势 。 
3.1.2 昆 仓 山 不 同 区 域 冰 湖 变化 特征 VA SITE 
界 将 昆仑 山 分 为 东西 两 部 分 ,分 别 讨论 各 区 域 冰 湖 
数量 和 面积 变化 情况 。 结 果 表 明 ,研究 区 冰 湖 数量 


一 兴 都 库 什 山 一 喜马拉雅 与 青藏 高 原 地 区 冰 湖 面 


与 面积 均 呈 西 多 东 少 的 分 布 特征 (图 4)。2020 年 ， 


y=12.464x+634.14 


80 m— 冰 湖 面积 /km? 一 4 一 冰 湖 数量 /个 R20.6242 + 1000 
---- 线性 ( 冰 湖 面积 ) 一 一 线性 ( 冰 湖 数量 ) 
+ 800 
60 
ge se 
gE y=1.1581x+27.188 nee 4 600 gg 
E 4 R=0.8659 — | 25 xd = 
x 2 898m - 400 3X 
20 
+ 200 
0 0 
2000 2005 2010 2015 2020 


年 份 


图 3 2000 一 2020 年 昆仑 山地 区 冰 湖 数量 与 面积 变化 统计 
Fig. 3 The number and area change statistics of glacial lakes in Kunlun Mountains from 2000 to 2020 
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(b) 东西 部 冰 湖 面积 
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图 4 2000 一 2020 年 昆仑 山地 区 东西 部 冰 湖 数量 与 面积 对 比 


Fig. 4 Comparison of the quantity and area of glacial lakes in eastern and western Kunlun Mountains from 2000 to 2020 
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西 昆仑 山 冰 湖 数量 占 总 量 的 73.60%, 冰 湖面 积 占 总 
面积 的 81.11%。 西 部 山区 海拔 高 ,冰川 宪 盖 面积 
阔 , 山 谷 间 常 有 积 雪 融 水 聚积 ,适宜 冰 湖 的 发 育 和 
扩张 。 而 东部 地 势 较为 平坦 ,多 发 育 有 以 降水 为 主 
要 补给 的 大 型 湖泊 。 从 冰 湖 数量 和 面积 变化 的 绝 
对 量 来 看 ,西部 冰 湖 的 数量 和 面积 增长 剧烈 ,成 为 
研究 区 冰 湖 数量 和 面积 增加 的 主要 成 分 。 近 20 a, 
西 昆仑 山 冰 湖 平均 以 每 年 11.6 个 的 速度 增长 ,东部 
年 均 冰 湖 增长 数 接近 于 0, 甚 至 有 大 量 年 份 出 现 负 
增长 现象 。 东 \ 西 昆仑 山 冰 湖面 积 变化 基本 一 致 ， 
西 昆仑 山 面积 年 均 增 长 率 为 426% , 东 昆 仑 山 面积 
年 均 增 长 率 为 3.78%。 总 体 而 言 , 西 昆 仑 山 冰 湖 所 
占 比重 大 ,增长 快 ,在 未 来 也 具有 很 大 的 扩张 潜力 。 
3.1.3 昆 仓 山 不 同 规模 冰 湖 变化 特征 ”由 于 人 研究 区 
范围 较 大 ,地 形 复杂 ,昆仑 山 冰 湖 规模 差异 大 ,面积 
最 小 的 约 0.01 km? (2000 年 ), 最 大 冰 湖 面积 可 达 
9.53 km? (2020 4E.) , UE 20 a 面积 平均 值 为 39.93 
km。 其 中 , 绝 大 多 数 冰 湖 为 面积 小 于 0.1 km 的 小 
型 冰 湖 。 

分 析 发 现 ,不 同 规模 的 冰 湖 数量 变化 趋势 具有 
一 致 性 , 均 与 冰 湖 总 数 的 变化 曲线 相 吻 合 ( 图 5$)。 
面积 小 于 0.1 km 的 冰 湖 数量 最 多 ,平均 占 冰 湖 总 数 
Ht WY 84.36% ,面积 超过 0.2 km? 的 冰 湖 数量 仅 占 
6.80% ,面积 越 大 的 冰 湖 ,数量 越 少 。 从 不 同 规模 面 
积 总 和 来 看 ,大 于 0.2 km 的 冰 湖 占 比 最 大 ,可 达 
55.17% ;其 次 为 面积 小 于 0.1 ker 的 冰 湖 ,平均 占 比 
为 26.03% ,0.1~0.2 km 的 冰 湖 面积 所 占 比 重 最 小 。 
面积 小 于 0.1 km 的 小 规模 冰 湖 数量 最 多 ,增长 速度 
最 快 ,数量 平均 增长 速度 为 每 年 13.45 个 ,演化 更 为 
激烈 。0.1~0.2 km 的 冰 湖 数量 近 20 a 表现 为 负增长 
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趋势 ,该 类 冰 湖 多 位 于 较 高 海拔 地 区 ,水 分 补给 受 
冰川 退缩 等 因素 的 影响 ,出 现 消退 甚至 消失 。 由 
此 ,昆仑 山 小 规模 冰 湖 演化 更 为 激烈 ,对 气候 及 其 
周围 环境 变化 敏感 ,与 Chen 等 ”认为 高 亚洲 地 区 小 
规模 冰 湖 (<0.04 km ) 增 长 最 为 明显 的 观点 具有 一 
致 性 ;大 规模 冰 湖 调 蓄 能 力 更 强 ,不 易 随 外 界 因素 
产生 变化 。 
3.1.4 昆仑 山 不 同 海拔 冰 湖 变化 特征 ”高 山地 区 的 
冰 湖 形成 多 与 冰 湖 周边 地 形 相关 ””。 为 了 进一步 分 
析 不 同 海拔 地 区 冰 湖 的 分 布 和 变化 ,本 研究 以 200 m 
间距 对 冰 湖 分 布 的 海拔 区 间 进 行 了 划分 (图 6)。 统 
计 结 果 显 示 , 近 20 a 中 国 昆 仑 山地 区 冰 湖 海拔 分 布 
的 下 限 为 3200~3400 m, 分 布 上 限 为 5800~6000 m, 略 
低 于 青藏 高 原 内 陆地 区 冰 湖 分 布 的 最 高 海拔 (2020 
青藏 高 原 冰 湖 分 布 最 高 海拔 为 62534m) ^, KE 
冰 湖 集中 分 布 于 4600~5600 mm ,数量 和 面积 分 别 占 
全 部 的 71.58% 和 70.51%。 这 是 由 于 人 研究 区 冰川 主 
要 分 布 于 海拔 4000~6000 m, 冰 川 面积 退缩 主要 和 集 
中 于 4500~5500 m, 超 过 此 范围 冰川 面积 儿 乎 无 变 
ft^. 昆仑 山 最 大 积 雪 量 位 于 海拔 5500 mm 左右 ,由 
于 积 雪 主要 来 源 于 高 海拔 地 区 的 固态 降水 ,可 推测 
人 研究 区 海拔 5500 m 处 存在 最 大 降水 带 近 。 充 足 的 
冰川 融 水 .固态 降水 等 补给 来 源 是 该 海拔 区 间 冰 湖 
分 布 较 高 的 重要 原因 。 
研究 区 冰 湖 数量 呈现 近似 正 态 分 布 ,在 4800~ 
5000 m 1t Fi] 4] fp Ec e ,平均 数量 为 137 个 , 约 为 总 数 
的 17.84%。 冰 湖面 积 主 要 分 布 于 5200~5400 m, P 
均 面 积 达 到 8.11 km , 约 为 总 面积 的 20.32%。 从 变 
化 趋势 上 看 , 冰 湖 数量 增加 集中 于 海拔 4800~5000 
m, 数 量 增加 了 46 个 , 约 为 总 增长 量 的 18.25% , 冰 湖 
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图 5 2000 一 2020 年 昆仑 山地 


区 不 同 规模 冰 湖 数量 与 面积 对 比 


Fig. 5 Comparison of the quantity and area of different scale glacial lakes in Kunlun Mountain area from 2000 to 2020 
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(a) 冰 湖 数量 
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图 6 2000 一 2020 年 昆仑 山地 区 不 同 海拔 冰 湖 数量 与 面积 统计 
Fig.6 Statistics on the number and area of glacial lakes at different elevations in Kunlun Mountains from 2000 to 2020 


面积 扩张 集中 于 海拔 5200~5400 m, 约 为 面积 增长 
总 量 的 20.32%。 

以 2020 年 为 例 , 对 昆仑 山 不 同 区 域 和 不 同 规模 
冰 湖 在 不 同 海拔 的 分 布 情况 进行 了 统计 分 析 ( 图 7)。 
对 于 不 同 区 域 冰 湖 而 言 , 东 昆仑 山 冰 湖 数量 分 布 峰 
值 位 于 海拔 4800~5$000 m, 冰 湖面 积分 布 峰值 位 于 
海拔 3600~3800 m, 在 此 区 间 内 ,东部 冰 湖 面积 为 西 
部 冰 湖 的 4.97 倍 。 造 成 以 上 特征 的 主要 原因 是 东 
昆仑 山北 坡 受 第 四 纪 以 来 地 质 断 陷 ,冰川 侵蚀 、 冰 
矿 阻 塞 等 的 影响 ,在 3600-3800 m 范 围 形成 了 多 
个 面积 较 大 的 冰 矿 阻塞 湖 ,导致 该 海拔 区 间 东 昆仑 
山 冰 湖面 积 异常 升 高 。 西 昆仑 山 冰 湖 数量 和 面积 
分 布 峰值 均 位 于 海拔 5200~5400 m, 向 两 端 呈 近 似 


对 于 不 同 规模 冰 湖 而 言 ,小 于 0.1 km! 冰 湖 数量 
和 面积 分 布 峰值 分 别 位 于 4800~5000 m 和 5000~ 
5200 m,0.1~0.2 km? 冰 湖 数量 和 面积 分 布 峰值 均 位 
于 5200~5400 m, 由 于 东 昆 仑 山大 面积 冰 研 阻塞 湖 
的 存在 ,大 于 0.2 km 冰 湖 面积 分 布 峰值 位 于 3600-~ 
3800 m。 小 规模 与 中 等 规模 冰 湖 易于 形成 和 相互 
转化 ,主要 由 地 质 历史 时 期 冰川 剥蚀 作用 产生 的 河 
槽 注 地 积 水 而 成 ,一 般 临近 冰川 /雪线 或 位 于 其 上 
方 的 较 高 海拔 区 域 。 大 规模 冰 湖 的 形成 对 于 地 形 
和 水 热 条 件 有 较 高 要 求 ,昆仑 山 雪 线 分 布 高 度 介 于 
4800-5300 m ^ ,印度 季风 和 太平 洋 季 风 越 过 青藏 
高 原 到 达 昆仑 山 时 影响 已 经 很 微弱 , 仅 可 形成 少量 
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降水 , 故 大 规模 冰 湖 多 分 布 于 相对 湿热 的 低 海拔 
区 域 。 

3.2 昆仑 山 冰 湖 分 布 变化 驱动 力 分 析 

32. 气候 变化 对 昆仑 山 冰 湖 变化 的 影响 近 20a 
昆仑 山地 区 气温 下 降 3.45% ,降水 减少 6.27% ,气温 
和 降水 数据 在 一 定 程 度 上 可 以 反映 该 地 区 气候 变 


2010 年 出 现 突 增 , 近 10 a 年 均 降 水 量 (159 mm) 较 前 
10 a(164 mm) 略 有 降低 。 由 降水 量变 化 趋势 图 可 直 
观看 出 , 东 昆 仑 山 降水 量 普遍 减 小 ,东南 部 少数 区 
域 呈 小 范围 不 显著 增加 趋势 ,而 昆仑 山东 部 的 冰 湖 
主要 集中 在 这 一 区 域 , 且 均 为 面积 较 大 冰 湖 。 西 昆 
仑 山大 范围 地 区 降水 量 保持 增加 ,在 78°*~80%E 降 水 


化 状况 (图 8)。 人 研究 区 气温 无 明显 变化 趋势 ,在 
2012 年 开始 上 升 。 气 温 变 化 趋势 图 直观 的 反映 了 
昆仑 山 气 温 的 时 空 演变 情况 (图 9) ,多数 区 域 气温 
无 显著 变化 ,东部 少数 区 域 呈 不 显著 增 温 ,稳定 的 
气温 可 为 冰 湖 的 形成 和 发 展 提供 条 件 。 西 昆仑 山 
部 分 区 域 气温 略 有 降低 ,可 能 导致 降雪 无 法 转化 为 
冰雪 融 水 ,从 而 限制 冰 湖 的 形成 ,使 该 地 部 分 冰 湖 
数量 和 面积 的 增长 率 有 所 下 降 。 

昆仑 山 降 水 量 呈 非 线 性 变化 ,2003 年 .2005 年 、 
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量 增 加 呈 极 显著 ,由 于 降水 的 大 幅 增 长 ,为 该 地 区 
冰川 发 育 提 供 了 充足 补给 ,间接 导致 冰 湖 面积 
张 ,面积 增长 率 逐 年 升 高 。 

3.2.2 冰川 变化 对 昆 仓 山 冰 湖 变化 的 影响 BTE 
仑 山地 区 气候 严寒 干旱 .降雨 稀少 RIRA ,冰川 
融 水 是 昆仑 山 冰 湖 补给 的 主要 形式 ”。 随 着 全 球 
气温 升 高 ,高 亚洲 地 区 冰川 快速 退缩 ” ,冰川 的 变 
化 成 为 影响 冰 湖 变化 的 重要 因素 之 一 。 本 人 研究 通 
it GEE 监督 分 类 提取 得 到 研究 区 冰川 ,并 与 中 国 西 
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图 7 2020 年 昆仑 山 不 同 区 域 与 不 同 规 模 冰 湖 海拔 分 布 统计 


Fig.7 Elevation distribution of different regions and glacial lakes of different sizes in Kunlun Mountain in 2020 
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图 8 2000 一 2020 年 昆仑 山 年 平均 气温 与 年 降水 量变 化 情况 


Fig. 8 Changes of annual mean temperature and annual precipitation in Kunlun Mountains from 2000 to 2020 
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北 冰 川 编目 数据 (2015 年 ) 进 行 对 比 , 确 保 结果 可 信 
后 统计 该 区 域 的 冰川 面积 变化 。 由 冰川 分 布 图 可 
知 (图 10), 昆 仑 山 冰 川 多 分 布 于 高 海拔 的 西 昆仑 山 
地 区 , 即 为 研究 区 冰 湖 数量 最 多 、 增 长 最 快 区 域 。 
冰川 面积 统计 结果 显示 , 近 20 a 昆仑 山 冰 川 面积 持 
续 退 缩 ,2000 年 昆仑 山 冰 川 面 积 为 12739.46 km? ,在 
2020 年 下 降 至 10045.51 km’, 共 减 少 21.15% o 
2000 一 2004 年 .2004 一 2009 年 冰川 面积 均 表 现 为 先 
升 高 后 缓慢 降低 的 趋势 。2010 一 2020 年 ,冰川 面积 
持续 降低 ,年 平均 减少 率 为 2.12%。 
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3.2.3 和 气温、 降水 和 冰川 面积 与 昆仑 山 冰 湖 的 相关 
性 “采用 皮尔 逊 相关 系数 法 探究 了 气温 .降水 和 冰 
川 面 积 与 总 冰 湖 及 不 同 区 域 .不 同 规模 .不同 海 拔 
冰 湖 之 间 的 相关 性 ( 表 1)。 结 果 表 明 , 冰 湖 与 气温 
和 降水 均 不 呈 显 车 相关 ,与 冰川 面积 多 表现 为 较 强 
的 相关 性 。 东 昆仑 山 距 离 冰 川 边缘 较 远 , 冰 湖 多 来 
源 于 冰川 模 谷 中 受 有 崩塌 ,滑坡 ,泥石流 等 地 质 灾害 
湾 堵 主 沟 道 而 形成 的 堰 塞 湖 ” 交 , 故 东 昆 仑 山 冰 湖 数 
量 与 冰川 面积 不 呈现 显著 相关 。 规 模 大 于 0.1 km? 
的 冰 湖 多 由 小 规模 冰 湖 间接 发 育 形成 , 故 与 冰川 面 
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图 9 Sen+Mann-Kendall 法 昆仑 山地 区 气温 .降水 变化 趋势 


Fig.9 Variation trend of temperature and precipitation in Kunlun Mountain area 
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图 10 昆仑 山地 区 2015 年 冰川 分 布 及 2000 一 2020 年 冰川 面积 统计 
Fig. 10 Glacier distribution map of Kunlun Mountain region in 2015 and statistical map of glacier area from 2000 to 2020 
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表 1 气温 \ 降 水 和 冰川 面积 与 昆仑 山 冰 湖 的 相关 关系 


Tab.1 Correlation of temperature, precipitation and glacier area with glacial lakes in Kunlun Mountains 


气温 与 气温 与 降水 与 降水 与 冰川 面积 与 冰川 面积 与 
冰 湖 数量 冰 湖 面积 冰 湖 数量 冰 湖 面积 冰 湖 数量 冰 湖 面积 

总 冰 湖 0.095 0.052 —0.085 —0.048 -0.573" -0.702" 

不 同 区 域 冰 湖 KEEL 0.178 0.079 -0.264 -0.021 0.031 -0.594" 
西 昆 仑 山 0.054 0.046 -0.013 -0.054 -0.687" -0.721" 

不 同 规模 冰 湖 <0.1 km? 0.114 0.001 -0.090 -0.021 -0.612" -0.729" 
0.1~0.2 km? -0.166 0.044 0.016 0.012 -0.401' -0.556” 

>0.2 km? 0.057 0.073 0.005 -0.070 -0.467° -0.673" 

不 同 海拔 冰 湖 3.2~3.4 km 0.115 0.039 -0.104 -0.137 -0.609" -0.697" 
3.4~3.6 km 0.107 0.036 -0.113 -0.039 -0.613" -0.708" 

3.6~3.8 km 0.073 0.051 -0.072 -0.048 -0.491' -0.702” 

3.8~4.0 km 0.196 0.001 -0.039 0.001 -0.509" -0.712" 

4.0-4.2 km 0.124 0.050 -0.135 -0.050 -0.603” -0.703” 

4.2-4.4 km 0.084 0.052 -0.076 -0.049 -0.579" -0.702” 

4.4-4.6 km 0.082 0.051 -0.098 -0.049 -0.575” -0.702” 

4.6~4.8 km 0.100 0.051 -0.088 -0.048 -0.562" -0.702” 

4.8-5.0 km 0.090 0.051 -0.085 -0.049 -0.573” -0.703” 

5.0~5.2 km 0.097 0.051 -0.079 -0.048 -0.572" -0.702” 

5.2~5.4 km 0.098 0.051 -0.086 -0.048 -0.566” -0.702” 

5.4~5.6 km 0.100 0.052 -0.078 -0.047 -0.584" -0.702” 

5.6~5.8 km 0.065 0.055 -0.062 -0.052 -0.549” -0.701" 

5.8~6.0 km -0.003 0.039 0.080 -0.069 -0.44T -0.683^ 


积 在 0.05 级 别 相关 性 显著 。 位 于 海拔 3600~3800 
m、3800~4000 m 及 5800~6000 m 的 冰 湖 数量 与 冰川 
面积 在 0.05 级 别 相关 性 显著 。 其 余 各 区 域 冰 湖 数 
量 和 面积 均 与 冰川 面积 在 0.01 水平 呈 显著 负 相 关 。 

昆仑 山 冰 湖 位 于 干旱 半 干 旱 高 山区 ,气温 和 降 
水 主要 通过 影响 冰川 面积 和 冰川 消融 速率 进而 影 
响 冰 湖 的 变化 ,其 与 冰 湖 变化 之 间 的 相关 性 均 不 显 
著 。 冰 川 融 水 作为 干旱 区 冰 湖 的 主要 补给 来 源 , 在 
气温 等 环境 因子 的 影响 下 ,冰川 加 速 融 化 导致 了 冰 
川 退缩 ,由 此 产生 的 大 量 冰 川 融 水 为 冰 湖 的 孕育 和 
扩张 提供 了 充足 的 补给 来 源 ,是 导致 新 疆 昆 仑 山 近 
20 a 冰 湖 数量 和 面积 迅速 增加 的 直接 原因 ,冰川 的 
变化 是 制约 冰 湖 发 展 的 重要 因素 。 

现 有 研究 表明 , 冰 湖 变化 是 气候 变化 、 冰 川 变 
化 湖 盆 环境 等 多 种 因素 综合 作用 的 结果 ,不同 地 
区 冰 湖 变化 的 主要 影响 因素 各 不 相同 。 喜 马 拉 雅 
地 区 和 天 山 山 脉 受 全 球 变 暖 的 影响 ,气温 升 高 导致 
区 域 冰川 快速 退缩 ,冰川 湖泊 不 断 增加 ”“”"。 念 青 
唐古拉 山 冰 湖 扩 张 速率 同时 受到 气温 和 降水 的 影 
Wap 。 青 藏 高 原 地 区 冰 湖 扩张 的 主要 驱动 因素 是 冰 


注 : 六 表示 在 0.01 级 别 ( 双 尾 ) 相 关 性 显著 ,* 表 示 在 0.05 级 别 ( 双 尾 ) 相 关 性 显著 。 


川 融 水 的 增加 ,同时 受降 水 量 和 流域 获 散 量 制约 ， 
与 本 研究 中 昆仑 山 冰 湖 变化 机 制 最 为 相似 。 


4 结论 


本 文 基于 GEE 平 台 进 行 监督 分 类 ,分 析 了 近 
20 a 来 昆仑 山地 区 冰 湖 面积 的 时 空 变化 特征 ,从 气 
温 、 降 水 和 冰川 面积 变化 等 方面 进行 驱动 力 分 析 ， 
结论 如 下 : 

(1) 2000 一 2020 年 ,昆仑 山地 区 冰 湖 数量 与 面 
积 总 体 表现 出 数量 增多 ,面积 增 大 的 趋势 ,但 增长 
幅度 有 所 下 降 。 截 至 2020 年 ,该 地 区 共有 冰 湖 数量 
894 个 ,总 面积 达 50.83 km。 冰 湖 数量 和 面积 增长 
量 分 别 为 39.25% 和 81.35% 。 研 究 区 冰 湖 数量 和 面 
积 均 表 现 出 西 多 东 少 的 特点 。 西 昆仑 山 冰 湖 数量 
平均 每 年 增长 11.6 个 ,和 东 昆 仑 山 年 均 冰 湖 增长 数 接 
近 于 0; 西 昆仑 山 面 积 年 均 增 长 率 为 4.26% , 东 昆 仑 
山 面积 年 均 增 长 率 为 3.78%。 

(2) 不 同 规模 的 冰 湖 对 气候 的 响应 不 同 ,面积 
小 于 0.1 km 的 冰 湖 对 气候 变化 更 敏感 ,增长 速度 最 


202308.00690v1 


ChinaXiv 


快 ,平均 每 年 增长 13.45 个。 面积 大 于 0.2 km 的 冰 
湖 占 比 最 大 ,其 次 为 小 于 0.1 km 的 冰 湖 , 占 比分 别 
为 55.17% 和 26.03%。 昆 仑 山 冰 湖 主要 集中 在 海拔 
4600~5600 m, 冰 湖 数量 增加 主要 在 海拔 4800~5000 
m, 面积 增加 集中 在 海拔 5200~5400 m, 

(3) 昆仑 山 冰 湖 数量 和 面积 与 气温 和 降水 均 无 
显著 相关 性 ,与 冰川 面积 呈 显 著 负 相关 ,表明 冰川 
面积 是 影响 冰 湖 变化 的 主要 因素 。 近 20 a 昆 仑 山 冰 
川 面积 退缩 21.15% ,冰川 退缩 带 来 的 积 雪 融 水 是 导 
致 冰 湖 数量 增加 ` 面 积 增 大 的 直接 原因 。 
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Abstract: Investigating the variable features of Kunlun glacier lakes is crucial for the development of the local 
ecological environment. The area and number of glacial lakes in the Kunlun Mountains have changed over the 
past 20 years, and this paper used supervised classification based on the Google Earth Engine remote sensing 
cloud platform to study the changes. It also examined the driving factors of temperature, precipitation, and glacier 
area. The findings indicate: (1) There were 39.25% more glacial lakes in the Kunlun Mountains in 2020 than 
there were in 2000, and the area expanded by 81.35%, with a distribution pattern of more lakes in the west and 
fewer lakes in the east. (2) The Kunlun Mountain glacial lakes are primarily found at an altitude of 4600-5600 m, 
and the number and area of glacial lakes account for 71.58% and 70.51% of the total, respectively. These glacial 
lakes have a smaller area than 0.1 km’ and are more sensitive to climate change. (3) The temperature and precipi- 
tation in the Kunlun Mountains declined by 3.45%, 6.27%, and 21.15% from 2000 to 2020, respectively, as did 
the glacier area. The primary cause of the expansion of glacial lakes is the meltwater produced by glacier melting. 
The study’ s findings may provide empirical justification for the preservation and use of water resources as well 
as for catastrophe warnings in dry regions. 
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